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In this paper we use the semiclassical Boltzmann equation to investigate the transport properties 

of Dirac ferniion on the surface of topological insulator with magnetic impurities. The results 

obtained show that there is also a minimal conductivity in this system as in graphene. We also 

argue the low temperature transport property, and find that there is no low temperature anomaly 

f^ ' known as Kondo effect when the temperature is T > 1Q~^ K. 
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CN ■ 1. INTRODUCTION 

Recently topological insulator has been theoretically predicted and experimentally observed in HgTe quantum 
q_) ' wells[l, 2], in Bii_xSbx alloys[3, 4], and in Bi2Se3 and Bi2Te3 bulk crystals [5-8]. The topological insulator is a new 
discovered novel material with gapped bulk phase and robust gapless surface(edge) state. And novel properties of 
topological insulator have been predicted, for instance, effective monopole and topological magnetoelectric effect [!)], 
superconductor proximity effect induced Majorana fermion states[10] etc. Sbi_xTex is the first material has been 
d • reported to be topological insulator[4], and Bi2Se3, Bi2Te3, Sb2Te3 have been predicted to be topological insulator[5] 
S [ with single Dirac cone on the surface. This material has potential applications in spintronics, topological quantum 
I I computation, eic[ll]. 

Recently, nanoribbons of Bi2Se3 have been fabricated [12], and transport properties of magnetic doped nanoribbons 

have been measured[13] with conductivity Kondo effect. However, the experiment can't distinct the conductivity 

Q [ anomaly from Kondo effect induced by bulk defect conductance. In this paper, we use the semi-classical Boltzmann 

equation to investigated the transport properties of topological insulator surface state with surface magnetic impurity 

doped. 

^ ' The paper is organized as follows: In Sec. II, we construct an effective interaction between topological insulator 

C^ i surface state and magnetic impurities, and give the corresponding Boltzmann equation with coherent distribution 
\l • function. In Sec. Ill, we analyze the incoherent case and give an exact solution for the Boltzmann equation. We also 
take into account the conductivity correction from electron-hole coherent terms in Sec. IV. In Sec. V, we consider 
the high-order correction of coupling constant, and give the renormalization group equation of the the interaction 
"i^ I between TISS and magnetic impurities. We make a conclusion in Sec. VI. 
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2. HAMILTONIAN AND BOLTZMANN EQUATION 



At first we take Sb2Te3 as an example to construct the Hamiltonian and make some assumptions: (i) The chemical 
5h I potential has been tune to seat at the Dirac point, which has been realized by D.Hsieh et aZ.[14] (ii) The interaction 
between topological insulator surface state(TISS) and magnetic impurities take the form of spin-spin interaction, (iii) 
The effective couplings are short-range rotational symmetric potential, which can be simulated by J^(r) = J^e~^/^/r, 
where exchange parameters J^ are estimated to the order O.leV ~ 0.5eV[15, IG], and the range R of interaction has 
been assume to about 13A[15]. (iv) Due to the effective RKKY interaction between magnetic impurities, all of the 
impurities have been ranged in the same direction[l-'')]. (v)The concentration of magnetic impurity is very small, so 
we can take the single impurity approximation. According to above, the Hamiltonian of this system can be described 
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by: 



Ho = hvF / d?rC\v)ia°'do.C{v), 



(1) 



F„t = J d'rJ^{r)St^CHr)a>'C{r), 



(2) 



where C(r),C^(r) are annihilation, creation operators of TISS respectively, cr" (a = .T,y) are the Pauli matrices in 
spin space of TISS, S^ [ji = x, y, z) is the spin operators of impurities and hvp — 3.7eV-A(wi? is Fermi velocity). The 
eigenstates of free Hamiltonian take: 
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with corresponding eigenvalues e±(k) = dzhvph and 9 is the azimuth angle aTctanky/kx- Under these basis, the 
scattering of TISS by magnetic impurity can be expressed as: 



/i,!/ — ±1 

Ct(ki)a(k2), 

-1/2 



(4) 



here J^ is the short-written of J^ (|ki — k2p + -tw) j and we assume J^ ~ Jy ~ J„. ^i, ^2 £^re the azimuth angle 

of ki and k2 . 

Electrical current of this system reads: 
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where /(k) is particle number distribution, and the velocity operator v is a 2 x 2 matrix in the space expand by 



V _ / Ga; cos 6* + Gj, sin 6* — i(Gj: sin6' — Gj, cos6') 
vp yi{exSmd — eyCos9) —ex cos 6 ~ By sin 6 



(6) 



here 9 is the azimuth angle of v. 

The Boltzmann equation, up to linear order of homogeneous electric field, of a similar system in graphene has been 

derived in reference [17]: 
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and the clastic collision term can be expressed by [18]: 
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where j^^ (Mj '^ = i) ^^e the four components of /(k), /g" are equilibrium distribution function without external 
electrical field, E is the homogeneous external electrical field, and v^j (i,j=l,2) are the four elements of velocity matrix. 



3. INCOHERENT CASE 



In equation (7) and (8), by setting the coherent terms / +(k),/+ (k) to be zero, we can get a more famihar 
Bohzmann transport equation of the form: 



eE • v„ 



j-^ (^- ) S{e,{k) e,.(k'))|(V^(k)|T|^^(k')>P (/^(k') - /,.(k)) , 
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where /i = ± for V''±(k) states, T is the transition matrix, E is the homogeneous external electric field. Under the 
first order of coupling constant T-Matrix can be expressed as: 



(u(ki)|f|w(k2)) = i [j,(e'(''i-''^) - 1)5^ + jy^'S- + J,e-'^'S^ 

and under the relaxation time approximation up to linear order of external homogeneous electric field, we have 

dfJk) 
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According to equation (9), (10) and (11), we can find an exact analytical solution for relaxation time: 

1 kR'^N 



r±(k) 2ffvFVTT4Wk^ 



JtS'S' 



1-3 



1 \ Vl+4i?2fc2_ (l + 2i?2fc2) 



2i?2fc2 



2i?2fc2 



+J'',{S{S + 1)-S'-S') 



^/l + 4J?2fc2 _ (1 _^ 2i?2fc2) 
^ 2^2p 



(12) 



here N is the concentration of magnetic impurities on topological insulator surface. In the short-range limit Rk <C 1, 
relaxation time takes: 
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4. COHERENT TERMS INCLUDED 

In this section, we take into account contributions of coherent terms in equation (7) and (8), the non-equilibrium 
part of the contribution function in these equations can be written as: 

/|+(k) = eEvF[cos{d-(l))Tn{k)+sm{9~(j))Xn{k)], 

fl_{k) = eEvF[-ism{9-(j))Ti2{k)~icos{0~q))Xi2{k)], 

/l+(k) = eEvF[ism{e~q})T2i{k)+icos{9^(j))X2i{k)], 

/l_(k) = e£;wF[-cos(6i-(/))T22(fc)-sin(6l-<?!))A22(fc)], (14) 

here is the azimuth angle arctan(fcj,/fc3;) and (j) determines the direction of electric field E = E(cos(/), sini/)). It can 
be verified that longitudinal current is proportional to the terms contain t^ (fc) and transverse current is proportional 
to the terms contain Xij(k). And we will show that the transverse current induced by magnetic impurities vanishes. 

We calculated all of the matrix elements K"g : 



K;^8 = T((^a(k)|r|V',(k')))((V/3(k)|f|V'.(k')))*, 



replaced non-equilibrium terms of f/ivO^) with equation (14), and took some tedious calculations from cquations(7) 
and (8), then wc found that Ty(fc) and \j{k) satisfy the following equations (in the short-range limit Rk <C 1): 



-2(a + C)Tii - a(rl2 + T21) + «C(Ai2 - A21) = 
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where a., £, and C arc three parameters: a = fH^J^S'S', ^ = fR'^[2J^S'S'+2J^{S{S+l)-S'S')], ( = -fR'^J^S'. 
(1 -f^ 2) means exchange indexes 1 and 2 to get another four equations. By solving Eq.(15), we find a solution to 
equations (7) and (8) in the short-range limit, which are: 
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here we define: 
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(16) 
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And the new relaxation time with coherent term correction is determined by 

hG/2 



-(fc).|e(fc)| = 
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(17) 



The conductivity of this system can be deduced from equation (5) straightforwardly, while the longitudinal conduc- 
tivity reads: 
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and the contribution to transverse current from terms in /ii(k) and /22(k) proportional to 1/H have been canceled 
to each other. 



HIGH ORDER CORRECTIONS 



Assuming the correction from coherent terms is small -p ^ 1, so the conductivity correction from high order terms 
of transition matrix may be taken into considered only in the incoherent case. In the short-range regime, Poorman's 



renormalization suggest that traditional RG cquation[] 0] is still correct in this system: 
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where D is the high energy truncation, with the magnitude of bulk energy gap, I is the renormalization rescaling 
factor: I = \og{D/T). Higher order correction of transition matrix elements have the form: 
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here we have considered an additional screened Coulomb potential in the interaction Hamiltonian, which has been 
induced by higher order interactions of TISS and magnetic impurities, so Jq ^ Jz,i\- It can be proved by following 
the Poorman's renormalization that the RG equation is dJo/dZ = 0. Finally the relaxation time correction can be 
calculated by the following equation: 
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And the corresponding conductivity correction is: 
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The integration over k' has been limited in a finite regime, because these conductivity correction is induced by a 
virtual state \Tpp{k')) while \ip^(k')) was propagating, and the uncertainty principle needs |ep(fc') — e^(k)\T{k') < h. 
Then according to equation (17), the boundary of ep{k') takes the form (.'i{k') = e^(fc)/(l ± ^). The conductivity 
correction can be rewritten as: 
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where k{2/G) (0 < 2/G < 1) is a nondimensional function of 2/G. So the high-order conductivity correction 6a 
proportional to temperature. In the low temperature limit, coupling constants flow into a strong coupled regime 
2/G — > 00, this semi-classical Boltzmann equation method is invalid because equation (17) is not consistent with 
uncertainty principle. For this reason, it is necessary to discuss the validity of our results obtained furthermore. 
According to equation(19), we estimate the Kondo temperature of this system and obtain Tk = -Dexp{ — b^jtj}^ 
where D take the value of bulk energy gap which is about O.leV for Sb2Te3, R is about 13A[15], and antiferromegnetic 
coupling is assumed to be isotropy with value 3.5eV. We get Kondo temperature is about 10~^K. So we obtain that 
surface magnetic dropped Sb2Te3 doesn't reveal low temperature conductivity anomaly when temperature is not very 
low and hope be observed in future. 



6. CONCLUSION 

111 this paper we investigate the transport properties of Dirac fermion on the surface of topological insulator with 
magnetic impurities. We find that there is also a minimal conductivity in the system studied as in graphene and there 
is no low temperature anomaly known as Kondo effect when the temperature T > 10"^-?^. 
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